16.451 Lecture 24 : The Nuclear Shell Model Dec. l, 2005 


If we choose the right potential f unction V N (r), then the wave f unction for the whole 
nucleus can be written as a product of the single particle wave functions for all A 
nucleons, or at least schematically: 

^ Nucleus ( F ) = n=l Vnlmin) 


oversimplification here... actually, it has to be written as an 
antisymmetrized product wavefunction since the nucleons are 
identical Fermions - the procedure is well-documented in 
advanced textbooks in any case! 

With total angular momentum given by: 

A 

J = Yjli> li = ~*i + S; , (s = j ) 

i= 1 


And parity: 



Always + for an even 
number of nucleons... 


What to use for V N (r)? - three candidate potential functions: 
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Square we) I | 

r< ft: V{r)~ - V g 
i l/fri=0 


Harmonicosciilator 

m = 


Saxon-Wood s 



Advantage: easy to write down 
Disadvantages: 

numerical solutions only 
edges unrealistically sharp 


Advantage: easy to write down and 
can be solved analytically 
Disadvantage: potential should 
not go to infinity, have to cut off 
the function at some finite r and 
adjust parameters to fit data. 


Advantage: same shape as measured 
charge density distributions of 
nuclei, smooth edge makes sense 

Disadvantage: 

numerical solution needed 
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Comparison: Harmonic Oscillator versus Woods-Saxon solutions: 3 


• since both potentials are spherically 
symmetric, the only difference is in 
the radial dependence of the wave 
functions 

• amazingly, when parameters are 
adjusted to make the average 
potential the same, as shown in the 
top panel, there is remarkably 
little difference in the radial 
probability densities for these 
two potential energy f unctions! 

• this being the case, the simplicity 
of the harmonic oscillator potential 
means that it is strongly pref erred 
as a model for nuclei 



Evidence that this works: 
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z 

electric charge density, measured via electron scattering: P$) = e Vi (0 1 2 

i = 1 


charge density difference 
between 205 TI and 206 Pb is 
proportional to the square 
of the wave f unction for 
the extra proton in 206 Pb, 
i.e. we can actually measure 
the square of the wave 
function for a single proton 
in a complex nucleus this 
way! 
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Solution: the “spin-orbit" potential 
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• Meyer and Jensen, 1949: enormous breakthrough at the time because it was the 
only explanation for the observed pattern of “magic numbers" and paved the way 
for a “periodic table" of nuclei ... and the Nobel prize in physics, 1963! 


• simple idea: 


(http://www.nobel.se/physics/laureates/1963/index.html 

- see Maria Goeppert-Meyer's Nobel Lecture link on this page) 


Vjv(r) => Vjv(r) + V so (r) 


£.s 


expectation value 


as in the calculation of magnetic moments, we can write: 


(i»s) - y(j 2_ ^ 2_s2 ) - ^L/(j+i)-^+i)-s(s+i)) 


• but there are only two ways the orbital and spin angular momentum can add for a 
single particle nucleon state: 


a) “stretched state" j - / + / : 

b) "jack-knife state" j = / - /: 





- s(/ + l) 


-(l + l) 
2 
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Spin-orbit force, continued: 
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• The energy shift due to the spin-orbit interaction is between states of the same 
/ but different j; 



\n,j = £- 1/2) 

\n,i) y 

5E 



|n, ] = 1+1/2) 


• the splitting is proportional to / and so it increases as the energy \r\crzaszs for the 

single particle solutions to V(r) 

• each state can accommodate (2j+l) neutrons or protons, each with different 

• empirically, the sign of the spin-orbit term for nuclei is opposite to that for atoms 

and the effect is much stronger in nuclei - the phenomenon has nothing to do with 
magnetism, which is the origin of this effect in atoms, but rather it ref lects a basic 
feature of the strong nuclear force. 

• with these features, the spin-orbit potential is the "missing link" required to 

correctly predict the observed sequence of magic numbers in nuclear physics © 
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Applications of the Shell Model 9 


Generic pattern of single particle states 
solved in a Woods-Saxon (rounded square 
well) potential model with appropriate 
spin-orbit interaction to reproduce the 
observed “magic number" pattern: 

State labels: J1 i j 

where n labels the order of occurrence of 
a given / value, and the state labels for 
orbital angular momentum are: 

0 1 2 3 4 5 6 

s p d f g h i 

Each state can hold (2j +1) neutrons 
and (2j+l) protons, corresponding to 
2(2j+l) distinct conf igurations of 
identical nucleons (m t , mj) to satisfy 
the Pauli exclusion principle 



Quantum numbers for a nucleus: 
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First of all, consider a "closed shell", which corresponds to a completely filled 
single-particle state, e.g. ls 1/2 , lp 3 / 2 , etc... containing (2j+l) protons or neutrons: 

The total angular momentum is: 


J = X ^ ’ M = X m i = 0 



all have the same each m value is different, 

j in a given shell ranging from -j to +j 


There is effectively only one conf iguration here, with total z-projection M = 0. 
Therefore, the total angular momentum of a closed shell must be J = 0! 


The total parity is: 


, = r 

i (-i)' = 

(-i)* 

2j +1 

= + : 



Since (2j+l) is always even, the parity of a closed shell is always positive. 
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Valence Nucleons: 
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For a closed shell + n nucleons, the angular momentum and parity is determined by 
the n "valence" nucleons, since the closed shell contributes J n - 0 + : 


n 



II 

M 

II 

(-1/ 

n _ 

z = 1 




The parity is uniquely determined, but there may be several different values of J 
that are consistent with angular momentum coupling rules. Residual interactions 
between the valence nucleons in principle determine which of the allowed J has the 
lowest energy - we can't predict this a priori but can learn from experiment. 


"Holes" — for a state that is almost full, it is simpler to consider angular momentum 
coupling for the missing nucleons than for the ones that are present: 


result for 
a closed shell 


n 

2 J+l 

- 1 1 


W 

ii 

to 

magnitudes of the two partial sums have 
to be the same, and M values opposite. 

Z =1 

z_(n+l) 



Total angular momentum for a collection of missing particles, i.e. holes, is the same 
as for that same collection of particles in a given shell model state. 


Magnetic Moments: 
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As in lecture 20, we can write: ^ = Sj J Mn = (/}»J 


(J+l) 


Predictions are only simple in the limit of one valence nucleon, and then we have: 

valence proton: 




= Us.p s + i ) 


Mn 


valence neutron: 


— 9 s,n S Mn 


Simplifying the results, we find: 


j = t + 1/2 9 = [aid- 1 / 2 ) + g s / 2 l/% 

i = -l _1/2 t Z ^ j Q +3/2 ) / (j +1 )-g s j/2(j+i)W 

with cj S p — + 5.58, 9sn ——3.83, — 1, — 0 
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Application with one valence nucleon: 17 0 
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_ 3/2 + 


5 


> 

3/2" 

i 


UJ . 



4 

^ 5/2" 


3 

__ 1/2 7 


2 



1 

— 1/2 + 


0 

5/2 + 4r " 



*809 


There are 8 protons and 9 neutrons, so we only 
need the low lying states in the shell model spectrum 
to understand the energy levels: 


lf7/2 
Id 3/2- 

ld5/2 

1 P 1/2 
1P3/2 

lSl/2 


(28) ^ 

f 

© 1 

1 

■2si/2 _ ; 

1 ® 1 

1 - 1 


: ® ' 

1 ' 


valence n 


} 


Ground state: 
full to here plus 
one neutron 


Ground state quantum numbers should be those 
of the valence neutron in the ld 5/2 state: 

J* = 5/2 + © 


Magnetic moment prediction: j = 1 + odd neutron 

^ = ^neutron = "1-91 H N 

measured value: -1.89 |i N excellent agreement! 


Excited states of 17 0 can be understood by promoting the valence neutron: 14 



First excited state: J n = } 



■ state 


ground state full to here 
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Excited states of 17 0: 
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Next excited state: J 71 = \ - 

explained by promoting a neutron from 
the filled 1 p 1/2 level to the ld 5/2 level 


lf7/2 - 
1^3/2- 


“ 251 / 2 ' 


■8 

.4 

‘2 



0 + pair 


neutron hole 


Excited states of 17 O 
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The 5/2' state is not so easy: to have negative parity, 
there must be an odd nucleon in a p state (or f state, 
but that is higher) 


lf7/2- 

ld3/2- 


-4 



But two neutrons are required to have different 
values of mj by the Pauli principle. Writing out 
the allowed conf igurations -> only J = 0, 2 , 4 
are allowed! (J = 2 will work here) 
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Two valence nucleons: (enough already!) 
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This problem is much more complicated! The inner "core" nucleons couple to J n = 0 + 
but in general there is more than one possibility for the angular momentum coupling 
of the valence particles. 

a) (pp) or (nn) case: Z and N are both even in these cases, so we know that 

the ground state conf iguration will be 0 + no matter what 
shell model state they are in. Excited states will have higher 
angular momentum, with possibilities restricted by the 
Pauli principle. 

b) (np) case: Z and N are both odd in this case. Only 6 examples in the whole 

nuclear chart!!! In isolation, (np) prefers to form a bound 
state - the deuteron - with J n = 1 + - 


U7/2 ■ 
Id 3/2. 


■ 251 / 2 - 


lDl /O Sr' — 


lPl/2- 
1P3/2 ■ 

1 S 1/2 " 




■8 

.4 

‘2 

"6 

- 2 ~ 

-4 


Example: 18 F, Z = 9, N = 9 

np pair: no restrictions on total J. 
possibilities are 0, 1, 2, 3, 4, 5 

Ground state of 18 F is 1 + i.e. deuteron 
quantum numbers! 

N.B. valence nucleons interact with each other, or 
all J values would be degenerate! 
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Odd neutron nuclei: 
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What is wrong? 

• the single particle model is too simple - nucleons interact with each other 

• conf igurations may be mixed, i.e. linear combinations of different shell model states 

• magnetic moments of bound nucleons may not be the same as those of free nucleons... 


Title page of a research monograph, Oxford, 1990: 


Modern Theories of 
Nuclear Moments 


B. CASTEL 

Qltrcrt\°. Vri Wttsiry. Cattiltiil 

AND 

I. S, TOWNER 

Chaik River Nuclear Cfiatk River, Canada 


From the preface: 
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Nuclear moments have had a long history which parallels the develop- 
ment of nuclear physics as a whole. In 1937, Schmidt proposed a standard 
for magnetic moments based on the assumption that the free-nuclcon 
magnetism persists in heavy nuclei. Then in the early 1950s attempts were 
made to explain deviations from the Schmidt estimates by invoking 
meson -exchange corrections or by modifying single -par tide wave func- 
tions so as to include many-body effects. The early 1970s saw the 
emergence of fundamental problems in trying to reconcile quadrupoie 
transitions and moments in light nuclei with shelimiodel predictions. 
Thus the study of nuclear moments could prove an excellent pedagogical 
tool to acquaint oneself with the complexities and elegance of some of 
the most current and powerful nuclear models. 
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Even the most sophisticated nuclear models are not completely successf ul! 21 



Fig. 2.1 Nuclear moments in the odd indium isotopes measured with laser 
spectroscopy techniques: (1) quadruple moments of / = 9/2* ground, states; (2) 
magnetic moments of / = 9/2* pound states; and (3) magnetic moments of 
/= 1/2’ isomeric states. The broken lines arc theoretical calculations based on 
core-particle coupling models (see Chapters 5 and 7) using effective coupling 
constants given in the figure. The Schmidt limit is shown as a dotted line. 


An interesting puzzle from beta decay: 
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last homework assignment: n C — ► n B + e + + v e log(Ft) =3.6 t i = 20.4 min 

14 C -> 14 N + e- + v log(Ft) = 9.0 t f = 5370 yr 
why so different? Both decays should have M if = 1 in the shell model (superallowed) 


U C 


14 C 




n n n n 


Q Q 


o Q 


neutrons 1 protons 
decay is fast, as expected 


ipi/ 
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o Q 




lPl/2 

- n - n - Q — ^ In 

i r3/2 


Q Q 


Q Q 


neutrons 1 protons 


ISt , 


decay is anomalously slow. - 
permits the use for 14 C dating! 


Explanation lies in interactions between nucleons that are not included in the shell model 
but required to fit experimental details. Even so, 14 C is not satisfactorily understood. 
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